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Experimental quantum information processing with superconducting circuits is rapidly advanc¬ 
ing, driven by innovation in two classes of devices, one involving planar micro-fabricated (2D) 
resonators, and the other involving machined three-dimensional (3D) cavities. We demonstrate that 
circuit quantum electrodynamics can be implemented in a multilayer superconducting structure that 
combines 2D and 3D advantages. We employ standard micro-fabrication techniques to pattern each 
layer, and rely on a vacuum gap between the layers to store the electromagnetic energy. Planar 
qubits are lithographically defined as an aperture in a conducting boundary of the resonators. We 
demonstrate the aperture concept by implementing an integrated, two cavity-modes, one transmon¬ 
qubit system. 


I. INTRODUCTION 

Circuit quantum electrodynamics (cQED) [TJ|2], based 
on the interactions of superconducting qubits with mi¬ 
crowave light, is currently emerging as one of the most 
promising experimental platforms for quantum informa¬ 
tion processing and quantum optics experiments Et 
[8]. In these superconducting circuits, Josephson junc¬ 
tions provide the non-linearity for qubits, while low-loss 
microwave resonators provide linear processing functions 
for quantum memories EHia, readout or entanglement 
buses [131 El 5 and filtering [T^ |T6] . 

Circuit QED has been developed in two platforms: 
fully planar (2D) circuits, which benefit from the geo¬ 
metric precision and parallel production of established 
micro-fabrication technologies, and 3D circuits involving 
conventionally machined cavities, but with superior co¬ 
herence times. Through improved design and material 
optimization, 2D qubits and resonators have significantly 
progressed, with internal quality factors (Qi) exceeding 
10^ (Refs. UMB). On the other hand, 3D resonators 
store a larger fraction of their electromagnetic energy in 
vacuum, making them less susceptible to material imper¬ 
fections, and can reach Q[ > 10^ (Refs. [TTl[^). Can a 
new cQED design take advantage of the benefits of both 
2D and 3D platforms? 

We propose to lithographically pattern qubits and res¬ 
onators in multiple planes separated by vacuum gaps 
used to store the electromagnetic energy. Thin-film alu¬ 
minum resonators built in this multilayer planar way re¬ 
cently demonstrated low-losses {Qi > 3 x 10^) at the 
single photon level [23]. One of the main challenges in 
the implementation of a multilayer approach to cQED is 
the design of qubit-resonator coupling between different 
layers of a structure. 

In the 2D and 3D platforms, coupling is achieved by 
inserting the qubit metallic structure onto the insulat¬ 
ing region of the resonator. In a perturbative descrip¬ 
tion of the coupling, the electric field of the resonator 
mode is aligned with the electric dipole of the qubit mode 
(Fig. la). However, in a multilayer architecture, this 
method would require fabricating qubits perpendicular 
to the lithographic planes. We propose a different de¬ 


sign strategy in which the qubit design layer coincides 
with one of the lithographic planes. It uses out-of-plane 
fields to couple the qubit—which we nickname aperture 
transmon—to the resonator mode (Fig. lb). 

To demonstrate the feasibility and advantages of our 
multilayer planar platform for cQED, we present in this 
letter the implementation and coherence properties of an 
integrated system composed of two standing modes cou¬ 
pled to a qubit, a now standard configuration for many 
basic cQED experiments [6] IH [24ll26] . Our implemen¬ 
tation is based on the two TEM modes of a supercon¬ 
ducting whispering gallery (WG) resonator introduced 
in Ref. [23]. One of the modes is over-coupled {Q = 10^) 
to a readout amplification chain, while the other is main¬ 
tained as high-Q as possible (Q = 2 x 10^). Both modes 
couple to a transmon qubit [27] with a Ti = 70 /as l ife- 
time. The Hamiltonian of this device (see Eq. Al) is 
similar to that of the 3D device in Ref. [6], and can be 
used for the implementation of cavity-based error correc¬ 
tion protocols uni. 


II. DEVICE AND METHODS 

Figure shows the multilayer chip-stack elements of 
the measured device. Two two sapphire chips served as 
substrates for each of the two Al patterned rings. We 
positioned the chips with the rings aligned inside an Al 
sample holder to establish the boundary conditions of 
the TEM modes of the WG resonator. Machined ledges 
in the sample holder maintained a 100 /am vacuum gap 
between the chips (details of the assembly can be found 
in the Appendix). The two orders of magnitude in as¬ 
pect ratio between the mode wavelength and the stack 
gap-spacing ensured tight confinement of the fields of the 
modes. In particular, the inductive participation ratio of 
the sample holder was found to be 10“^ or smaller for 
each of the modes, as computed with an HFSS [28] finite 
element model, assuming a London penetration depth of 
50 nm. 

The transmon qubit was directly patterned in the thin- 
film of the ring in layer 1, as shown in the inset of Fig. 2. 
The qubit consists of a 0.05 x 0.5 mm island inside a 
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FIG. 1. Qubit-resonator coupling in different cQED ap¬ 
proaches. a) In-plane coupling in 2D. The electric field lines 
of the resonator (blue) are aligned with the dipole moment 
of the qubit (red), both of which are in the plane of qubit 
fabrication, b) Out-of-plane coupling in a multilayer planar 
device. The resonator is now represented as a section of a 
multilayer whispering gallery mode resonator consisting 
of two superconducting thin-hlm rings deposited on different 
sapphire substrates that are separated by an electrically-thin 
vacuum gap. The qubit is defined by an aperture carved di¬ 
rectly from the conducting boundary of the resonator. The 
orange and blue arrows represent the resonator surface cur¬ 
rent density and electric field lines, respectively. 

0.23 X 1 mm aperture in the ring, connected through a 
Josephson junction with Ej/h = 12 GHz. The qubit 
structure perturbs the resonator mode frequencies only 
at the percent level. The island nominally shares 60 fF 
of capacitance with its own ring and 5 fF with the oppo¬ 
site ring. The junction capacitance and these geometric 
capacitances define the nominal qubit charging energy 
Ec/h = 275 MHz, frequency ujq^l27T = 4.85 GHz, anhar- 
monicity a = 320 MHz [27]. The rings and the qubit 
are fabricated simultaneously in a single electron-beam 
lithography step using a double-angle, bridge-free tech¬ 
nique [2^ [30] . 

The spatial mode orthogonality of the two WG modes 
allowed us to implement the long lived storage (D-^) and 
over-coupled readout {D^^) modes of a quantum regis¬ 
ter within the same physical structure. In the follow¬ 
ing, we refer to these two modes simply as “storage” 
and “readout,” with nominal coupling > 10^ and 
= 1.8 X 10^ (see Appendix [a|) , respectively. 

The sample holder was thermally anchored to the 
15 mK stage of a dilution unit. We used the standard 
cQED measurement setup (See Appendix]^ with the ad¬ 
dition of a phase preserving, quantum limited, Josephson 
parametric amplifier [31]. 

III. PLANAR TO NON-PLANAR COUPLING 

For dispersive coupling between a transmon qubit (q) 
and a resonator mode (r), the strength of the cross-Kerr 
Xqr depends primarily on the detuning A = aper¬ 

ture geometry, and the resonator fields at its position. In 
order to quantify the aperture coupling independently of 


FIG. 2. Photograph of chip-stack elements. Thin-film A1 
rings are patterned in a single e-beam lithography step along 
with the Josephson junction on a sapphire substrate. The 
boxed region shows a magnified optical micrograph of the em¬ 
bedded aperture qubit (false-colored in red) and the location 
of the Josephson junction (red cross). The axis of symme¬ 
try, represented over substrate layer 2, defines the angular 
position 0 around the ring. 

the potentially tuna ble A, we d efine the effective cou¬ 
pling rate ^qr = A^Xqr(A)/£’c El], which is approxi¬ 
mately independent of A and Ej in the transmon limit 
EjjEc ^ 1 and for weak interaction ^ <C A. 

In 2D and 3D, the coupling strength g can be under¬ 
stood as arising from an interaction between the bare 
electric field of the resonator and the electric-dipole-like 
charge distribution of the qubit. Here, in our multilayer 
structure where the qubit is patterned in an aperture in 
one of the layers, the coupling mechanism is more in¬ 
volved. Figure and show the charge and current 
distribution of the qubit mode, respectively. The cou¬ 
pling is determined by the overlap between these distri¬ 
butions and those of the resonator mode. 

The interplay of the capacitive (charge overlap) and 
inductive (current overlap) coupling is shown in Fig. W’ 
where we plot the dependence of g on the qubit position 0 
(see Appendix H for the simulation procedure based on 
black-box circuit quantization [32]). Varying the position 
0 of the transmon varies its coupling g to the resonator, 
independently of its frequency and anharmonicity. For a 
given position 6>, g can be further adjusted by changing 
the dimensions of the aperture. 

IV. EXPERIMENTAL RESULTS 

Microwave spectroscopy revealed the transmon qubit, 
storage (7^^), and readout modes at 4.890 GHz, 

7.070 GHz, and 7.267 GHz, respectively, in 1% agreement 
with the HFSS numerical simulations of the sample [33] . 
From qubit spectroscopy, we observed a 310 MHz trans¬ 
mon anharmonicity, in 5% agreement with the nominal 
qubit charging energy EH- 

Figure]^ shows the qubit free decay with an exponen¬ 
tial time constant Ti = 70 //s. We measured = 8 gs 
and T 2 = 20 //s. The dephasing noise was measured by 
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Electric field of resonator modes Aperture qubit mode Qubit-cavity coupling 



FIG. 3. Electric and magnetic contributions to qubit-resonator coupling, as calculated using HESS, a) Electric field amplitude 
|E I in the resonator WG modes as a function of 0, for one photon of energy. The angular dependence of the surface-current- 
density js follows that of E, but is shifted by 90°. b) Surface-charge-density amplitude a shown in color scale with overlayed 
white signs to indicate the relative charge polarity. The charges in the island and corresponding image charges in the opposite 
layer below determine the electric contribution to the qubit-resonator coupling, c) Qubit mode surface-current amplitude js 
shown in color scale with overlayed white arrows to represent the direction of the flow. The narrow rails on each side of the 
aperture are equivalent to a shared inductance between the qubit and resonator and determine the magnetic contribution to the 
qubit-resonator coupling, d) Qubit-resonator coupling rate g of the aperture qubit to the WG resonator modes as a function of 
the qubit position around the ring 6. The right vertical axis also shows the equivalent cross-Kerr y for a A = 1 GHz detuning. 
The coupling g is the algebraic sum of the electric and magnetic contributions, which interfere constructively or destructively 
as a function of 6. Eor the readout (storage), maximal constructive (destructive) interference occurs at about 45°, while near 
^ = 0,180° the coupling is capacitive (inductive). 


a CPMG technique (see Fig. A2). The contributions to 
the dephasing noise are not currently understood, they 
could be the result of photon shot noise [34], mechanical 
vibrations and/or offset charge drifts [27]. The readout 
linewidth, = 0.35 MHz, and qubit dispersive shift, 

Xqr/27r = 0.30 MHz, agreed to 10% with simulations. 


From spectroscopy of the storage at photon number 
n ~ 10^, we measure a linewidth Aa;s/27r = 4 kHz 
which includes self-Kerr broadening in addition to de¬ 
phasing. From this linewidth, we infer an approxi¬ 
mate lower bound on the storage coherence time: ^ 

2/Aco’s = 80 /is. 


To measure the storage lifetime Tf in the single-photon 
regime (see Fig. [^), we used the photon-number parity 
protocol introduced in Ref. [6l[25l[35], which, in our case, 
is more sensitive than a direct amplitude decay measure¬ 
ment. A 1 /is Gaussian pulse first displaced the stor¬ 
age cavity to a coherent state with n = 2.5 photons, a 
state with essentially zero parity Pg = (exp(—iTra’^'a)) = 
5 X 10“^, where a is the storage lowering operator. After 
a variable delay, a Ramsey-like sequence with a fixed time 
delay of n/xqs mapped the parity of the storage photon 
number to the qubit state. From the parity measurement, 
we extracted a low-photon number Tf = 45 /is. Measure¬ 
ments at higher photon numbers (up to n = 200) showed 
no power dependence of Tf. The calibrations needed for 
this protocol are described in Appendix [Pj 


V. DISCUSSION 


We summarize the two cavity modes, one-qubit de¬ 
vice interaction strengths and lifetimes in Table 1. The 
measured frequencies and coupling energies of the mul¬ 
tilayer device agree at the percent and ten-percent level, 
respectively, with design values from numerical simula¬ 
tions. The discrepancy can be explained by machining 
tolerances (25 /im) of the gap spacing and chip align¬ 
ment in the sample holder, and could be improved by 
using micro-machined separators to support the struc¬ 
ture [36H38] . 

The measured coherence times are on par with those of 
two-cavity, one-qubit devices using 3D rectangular cavi¬ 
ties liiHiEaEaiMiian]- The qubit and storage lifetimes 
were not limited by their input-output (I/O) conpling [AT| 
These lifetimes could be extended by design optimiza¬ 
tion, as well as material advances demonstrated in the 
2D and 3D architectures miiiiiiiHiisiiiiniiii]. Spurious 
fringing fields in the substrates and environment would 
be reduced by decreasing the gap spacing and improv¬ 
ing the chip stack alignment. In the present device, the 
100 fim gap captures ^90% of the cavity energy and 15% 
of the qubit energy in the vacuum. A decrease in the gap 
by a factor of 10 would, for both modes, decrease the 
bulk dielectric participation down to the percent level, 
a gain of more than one order of magnitude over planar 
and 3D qubits, provided that we would not be limited by 
the surface quality of the superconducting film. 

We demonstrated at least three orders of magnitude 
separation in I/O coupling Q between the storage and 
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FIG. 4. a) Aperture transmon energy relaxation. A single 
exponential fits the data with a Ti = 70 //s. The popula¬ 
tion inversion is defined as (Pe — Pe^)/(Pg^ — Pe^), where Pj^ 
and Pg^ are the thermal populations, b) Storage parity re¬ 
laxation from which we infer an energy relaxation lifetime of 
Ti = 45/iS. Inset: pulse sequence for the measurement, see 
Appendix 


Mode 

Qubit 

Storage 

Readout 

Frequency (GHz) 

4.890 

7.070 

7.267 

Tl (ytiS) 

70 

45 

0.42 

T2 I Ti ins) 

8* 1 20 

>80 

- 

ajlTT (MHz) 

310 

1 X lO-^ 

1.5 X 10“^ 

Xq/27r (MHz) 

- 

0.25 

0.30 


TABLE I. Main parameters of sample. The cross-Kerr in¬ 
teraction with the qubit mode is denoted Xq, while a de¬ 
notes anharmonicity. All parameters are measured except 
the storage/readout anharmonicity, which are calculated from 
<^s,r = Xq s,r/4<aq [32]. Thc symbol * indicates a Gaussian de¬ 
cay. 


readout, which not only share the same physical foot¬ 
print, but differ by only 200 MHz in frequency. This 
type of spatial mode control is advantageous when deal¬ 
ing with cross-talk and frequency crowding in devices 
with increased complexity. 

The measured device is a suitable candidate for a 
quantum register [TUI 112, with storage coherence time 
Ti = 80 /is exceeding that of the qubit by an order of 
magnitude. The storage could provide a large Hilbert 
space to encode quantum information, while the long- 
lived qubit serves as a conditional, non-linear control over 
the cavity space with a low bit-flip error rate. 

The qubit-resonator coupling geometry presented in 
Fig. lb relies on the use of an aperture in one of the 
metal layers. Radiation fields from an aperture usually 
constitute spurious loss and cross-talk mechanisms, but 
in our case, this effect is mitigated by the proximity of 
the opposite superconducting layer. In fact, our work 
demonstrates that we can put these fields to a good use: 
mediating the coupling between the planar qubit and 
multilayer resonator. This approach can be extended to 


provide low cross-talk inter-layer connections for devices 
with more than two layers, such as the architecture pro¬ 
posed in Ref. [38] . 


VI. CONCLUSIONS AND PERSPECTIVES 

We have implemented a multilayer superconducting 
device for quantum information processing that com¬ 
bines the benefits inherent to the precise geometry con¬ 
trol of 2D micro-fabrication with those of the coherence 
in 3D qubits and resonators. In particular, the qubit- 
resonator mode couplings can be precisely adjusted. We 
believe that the quality of the measured coherence in the 
present work resulted from the confinement of electric 
fields within the vacuum gap separating lithographically 
defined layers. The design principles illustrated by our 
work can be extended to devices with more than two 
layers, each layer corresponding to a specific function: 
qubits, control lines, resonators, amplifiers, etc. In par¬ 
ticular, the aperture based coupling method introduced 
here can be generalized to inter-layer coupling in such 
multilayer devices. Furthermore, the TEM mode struc¬ 
ture and the separation of layers provides a favorable ge¬ 
ometry for hybrid systems, such as spin-ensembles with 
cavities [43] [44| , spin qubits with magnetic contacts [45] , 
or nano-wire qubits [46] [47]. 
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Appendix A: Design details 


System Hamiltonian. If we expand the Josephson 
junction cosine potential to fourth order, apply the ro¬ 
tating wave approximation, and limit the Hilbert space 
of the transmon mode to the first two levels [32], then 
the effective device Hamiltonian is: 

1-LlTi =^(1 + ^ z ) + + ujrl)^h 

- +(^2.){Xq_sa^a + Xqibh) - Xstd'^ ab'^b, 

(Al) 


where a and b are storage and readout bosonic opera¬ 
tors, respectively, and ctz is a qubit Pauli operator. The 
storage-readout cross-Kerr is Xsr ~ XqsXqr/o^q, where o^q 
is the transmon anharmonicity. 

I/O coupling. Figure Air shows a simulated field 
profile for the DH mode, where the maximum currents 
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flow parallel to the axis of symmetry of the rings [23] . 
The orthogonal mode corresponds to exchanging the 
field and current maxima and minima. 

As illustrated in the inset of Fig. |A1| 3, two non¬ 
magnetic pins penetrated the sample holder lid to cou¬ 
ple capacitively with the maximal charge densities of 
the readout above the thinnest (0 = 0) and thickest 
{0 = 180°) parts of the rings. Owing to the selective cou¬ 
pling due to the spatial mode orthogonality, the nominal 
readout coupling was 1.8 x 10^, while the nominal 
storage mode coupling exceeded 10^. 


(a) Differential parallel (D") mode 

E 




I/O pin depth (mm) 


FIG. Al. a) HFSS calculation of the surface currents and 
electric fields, which are contained within the 100 /im vac¬ 
uum gap separating the two layers, for the TEM WG 
mode. The orthogonal mode closely resembles the 
up to a 90° rotation. Depending on the configuration of the 
input/output (I/O) coupling pins (see c) these modes can be 
used either as a qubit readout mode (low coupling-quality- 
factor Qc) or storage for quantum states (high Qc)- b) Sim¬ 
ulated Qc as a function of the coupling coax pin depth inside 
the sample holder at ^ = 0°. The inset shows a not-to-scale 
cross-sectional representation of the chip stack in a sample 
holder (black). We can selectively couple to (Q^ = 10^) 
while remaining uncoupled from (Q^ > 10^)? as indicated 
by the vertical gray line, which corresponds to the nominal 
parameters of the measured device. 


in a single lithography step on a PMMA/MAA resist bi¬ 
layer. We then performed a double angle Al evaporation, 
20 and 30 nm thick, in a Plassys UMS300 at a pressure of 
5 X 10“^ Torr. Between these two depositions, an AlOx 
barrier was formed by thermal oxidation for 6 minutes 
in a static 100 Torr environment of 85% argon and 15% 
oxygen. Chips were diced to 15.5 x 25.4 mm. 

Qubit design details. The qubit island inside is con¬ 
nected to the ring by a 1 /im thin wire and a 130 x 700 nm 
Josephson junction with Ej/h = 12 GHz. The value 
E:sIEq = 44 yields a maximum offset-charge dispersion 
of 30 kHz. 

Sample holder. The chips are placed inside the bot¬ 
tom piece on ledges that are machined in the Al walls. 
The sample holder top piece has four legs which use in¬ 
dium to secure the chips against the sample holder bot¬ 
tom. 

Measurement setup. An aluminium and permalloy 
shield protected the sample from stray magnetic fields. 
The SMA input lines had thermalized cryogenic atten¬ 
uators (20,10,30 dB) at the 4 K, 0.1 K and 15 mK 
stages of a Cryoconcept DR-JT-S-200-10 dilution re¬ 
frigerator, respectively. The sample holder output con¬ 
nected to a Josephson parametric converter (JPC) ampli¬ 
fier [31] through two Pamtech 4-8 GHz cryogenic circu¬ 
lators and superconducting NbTi-NbTi coax cables. The 
JPC served as a phase preserving amplifier which oper¬ 
ated near the quantum limit with a gain of 21 dB over 
a bandwidth of 5.6 MHz. Two 4-8 GHz circulators to¬ 
gether with two low-pass filters—a 12 GHz K&L multi¬ 
section lowpass and a box-type Eccosorb CR-110 filter— 
serve to isolate the JPC from the following Low Noise 
Factory HEMT with 40 dB of gain. We found an 8 dB 
noise rise figure for the amplification chain, indicating 
that the observed noise at room temperature is ^90% 
percent amplified quantum fluctuations, though the to¬ 
tal quantum efficiency of the measurements is lower due 
to losses. At room temperature, a 30 dB Miteq amplifier 
further amplifies the signal and feeds it into a standard 
heterodyne microwave interferometer operating at an in¬ 
termediate frequency of 50 MHz. An analog to digital 
converter records the mixed-down output signal together 
with a mixed-down reference of the input signal. Com¬ 
bining the output and reference signals accounts for any 
phase drift in the readout control generator. 

At room temperature, we used a Tektronix 5014C ar¬ 
bitrary waveform generator, an Agilent E8257D vector 
generator, and several Vaunix Lab Brick generators to 
generate the qubit and cavity tones. 


Appendix C: Details of qubit and storage mode 
Appendix B: Experimental details coherence measurements. 


Fabrication. We micro-fabricated both layers on the 
same 430 /im thick, double-side-polished, c-plane sap¬ 
phire wafer. Using a 100 kV Vistec electron beam pat¬ 
tern generator, we defined the WG resonator and qubit 


The qubit parameters were in the convenient regime for 
continuous state monitoring, where the dispersive shift 
and output coupling rate are nearly equal. From quan¬ 
tum jump measurements, not presented here, we inferred 
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a qubit excited state population below 4%. 



N (Number of echo pulses) 

FIG. A2. We measured = 8 /is (inset) and T 2 = 20 /xs. 
By dynamically decoupling the qubit from low-frequency 
noise we observed an order of magnitude improvement in the 
coherence time of the aperture transmon exceeding Ti. The 
dashed line is a guide for the eye. 


is more sensitive than a direct amplitude decay measure¬ 
ment. The dispersive readout is based on the cross-Kerr 
frequency shift of the readout mode due to the storage 
photon occupation. We apply a 500 /xs coherent pulse 
to excite the storage mode, followed by a 250 /xs, low- 
power tone to probe the readout frequency. The storage 
photon population decayed exponentially with a lifetime 
Tf = 45 /XS. The resolution of this dispersive measure¬ 
ment is too low to access the single-photon regime, be¬ 
cause of the small ratio (10“^) between the cross-Kerr 
and the readout linewidth. 

I/O coupling. Over several cooldowns, we progres¬ 
sively decreased the I/O coupling by shortening the cou¬ 
pling pin lengths. However, while the over-coupled read¬ 
out lifetime increased by a factor of two to the value in 
Table 1, the qubit and storage mode lifetimes, as well as 
the frequencies and non-linear coupling strengths, did not 
change measurably. From this, we place a lower bound 
on the measured storage coupling > 10^. 


Qubit coherence. The inset of Fig. |A2| shows the 
Ramsey coherence signal of the qubit which decays with a 
Gaussian envelope and a time constant = 8 /xs. Since 
this Gaussian envelope is indicative of low-frequency 
noise, we used dynamical decoupling techniques to access 
the intrinsic qubit coherence. A Carr-Purcell-Meiboom- 
Gill (CPMG) protocol, following the approach and pulse- 
train calibrations of Ref. [48], shifted the maximum of 
the longitudinal noise-susceptibility of the qubit to higher 
frequencies. Figure [A^ shows the increase of the dynam¬ 
ically decoupled coherence time T 2 as a function of the 
number of CPMG pulses. The improvement of T 2 be¬ 
yond Ti confirms the dominance of low-frequency noise. 
The Gaussian envelope, indicative of low-frequency noise, 
could be the result of mechanical vibrations and/or offset 
charge drifts. 


Appendix D: Photon number parity calibration with 
qubit-state revivals 


Qubit state revival 



Ramsey delay (|js) 



FIG. A3. Using the dispersive cross-Kerr interaction with 
the readout mode Xsr/27r = 0.25 kHz (pulse sequence in in¬ 
set), we measured a storage lifetime Tf = 45 /xs for h > 10 
photons. 


Storage mode lifetime. To measure the storage Tf 
at the 10-200 photon level, we u sed a dispersive readout 
of the storage ring-down (see Fig. A3), which, in our case. 


FIG. A4. Galibration of the photon number parity measure¬ 
ment in the storage was achieved with a qubit-state revival 
experiment. For small storage mode displacements fig ~ 0 
(blue), the decay is dominated by intrinsic qubit decoher¬ 
ence. For increasing displacements, up to fXg = 3 (red), the 
apparent increase in decoherence is due to the large qubit- 
cavity interaction rate, and we observe qubit-state revivals at 
integer multiples of 2/p = 27r/xqs. 

Experimental method of parity measurement 
and calibration. To measure the storage photon- 
number parity presented in Fig. 4b, we perform a pro¬ 
tocol introduced in Refs. 123 ES]. After displacing the 
storage using a coherent drive, we apply a 7r/2 pulse, 
which activates the qubit-storage cross-Kerr interaction 
and fully entangles the storage parity with the qubit at 
a time tp after the qubit pulse. A second 7r/2 pulse maps 
the parity to the expectation value of the qubit op¬ 
erator, which is read out projectively. This is the mea¬ 
surement sequence used for the parity measurement of 
Fig. 4b. 
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This protocol requires calibration of the initial dis¬ 
placement photon-number n and the parity mapping de¬ 
lay time tp = Tr/xqg. To perform the calibration, we 
displaced the WG storage by a short, coherent drive, 
and then performed a standard qubit T 2 Ramsey experi¬ 
ment, as shown in Fig. |A4| Sharp coherence peaks stro- 
boscopically reappear at integer multiples of the cross- 
Kerr interaction period 27r/xqs, indicating the value of 
Xqs/27T_ = 0.25 MHz. From a global fit to the theory 
(Eq. |D1| ) over all displacement amplitudes, we calibrate 
the corresponding storage photon numbers n. 

Calibration theory. During the measurement, the 
readout mode remains unpopulated, and we can ignore 
its contribution to the system Hamiltonian from Eq. |A1| 
In the rotating frame of the storage and qubit, the system 
Hamiltonian takes the form: 

H/h = -Xqsa'a\e) (e| . 

Eor a system starting in the ground state, the calibration 
Ramsey signal of Eig. |A4| obeys the following form as a 
function of time t: 

(cT,) = 

(Dl) 


where n is the average photon number in the storage 
mode, A is the pulse detuning from the qubit frequency, 
and I/T 2 is the incoherent dephasing rate. 


Appendix E: Simulation of qubit—cavity coupling 

We numerically simulated the qubit design shown in 
Eig. 2a for various qubit-position angles 0 using HESS, 
and for each simulation, we extracted the effective cou¬ 
pling rate g. We treated the Josephson junction as a 
lumped, linear inductor in each HESS eigenmode simu¬ 
lation and found the linearized mode frequencies to con¬ 
struct the linearized system Hamiltonian [32]. To treat 
the perturbing effect of the non-linear Josephson terms 
in the full Hamiltonian, we first calculated their mag¬ 
nitude using the energy participation-ratio method [49] . 
which is based on the fields already found in the eigen¬ 
mode simulation. Second, we numerically diagonalized 
the full Hamiltonian to find the energy spectrum of the 
system. Erom the spectrum, we extracted the frequencies 
and Kerr coefficients of the transmon and the resonator 
modes; from these, we calculated the coupling rate g. 
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